Priedel sum rule for an interacting multi-orbital quantum dot 
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A generalized Friedel sum rule is derived for a quantum dot with internal orbital and spin degrees 
of freedom. The result is valid when all many-body correlations are taken into account and it links 
the phase shift of the scattered electron to the displacement of its spectral density into the dot. 
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The Friedel sum rule (FSR) is one of the few exact 
results of solid state physics m , with a vaste range of ap- 
plications in the fields of scattering theory, magnetic and 
non magnetic impurities in metals 1], Kondo effect 0, 
and, very recently, coherent transport through quantum 
dots (QDs) II and molecules 0. This powerful relation 
connects the charge density displaced by an impu- 
rity or a nano-object, acting as a scattering center in a 
conductor, with the variation of the phase shift 6p of the 
scattered wave with respect to the energy huj: 
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The Friedel phase of Eq. (Q, which holds for single- 
channel elastic scattering only, is linked to the eigenvalue 
of the scattering matrix S through the identity 5' = e^'*'' . 
In the many-channel case, appearing on the left hand 
side of is replaced by TrlnS'/27ri p\. 

Recently Lee % and Taniguchi and Biittiker JJ showed 
the relevance of FSR in measurements of the transmission 
phase acquired by an electron passing through a QD em- 
bedded in the arm of an Aharonov-Bohm interferometer 
. These experiments |^ |^ allow for both directly mea- 
suring the phase shift and arbitrary controlling the Fermi 
energy (or, equivalently, the plunger gate voltage of the 
QD), namely the two quantities appearing on both sides 
of . Even if the identification of 5y with the transmis- 
sion phase is unjustified in generic situations 6, 7], nev- 
ertheless Aharonov-Bohm interferometry paves the way 
to the direct experimental test of FSR for a QD whose in- 
ternal structure, charge, spin, and correlation can be ex- 
ternally controlled. Indeed, the "non-universal" behavior 
of the QD transmission phase in the regime of very few 
electrons {N < 10) suggests that the QD orbital and spin 
degrees of freedom may play a major role ;9]. 

In real QDs used in interference experiments single- 
particle levels have a small energy separation (~ 0.5 
meV), if compared to characteristic Coulomb energies 
(~ 1 — 3 meV) j^, and therefore many of them should 
be included in any reliable model for electron correlation 
[lol im . The FSR is generally believed to hold even in 
the presence of electron-electron interaction . This was 
rigorously demonstrated only in two cases: (i) electrons 
in the metal form a Fermi liquid and the impurity has no 
internal degrees of freedom pj (ii) the interaction is lim- 



ited to a localized impurity orbital (Anderson model 
and its specific extensions to open atomic shells 0, F 
This Letter shows that the FSR must be reformulated 
in the experimentally relevant case of a multi-level in- 
teracting dot. A generalized statement holds, with the 
spectral density Af{uj) of the scattered electron accumu- 
lated in the QD replacing the non-interacting density of 
states g{u!)/e appearing in Q. We discuss the relevance 
of this result for the qualitative understanding of a few 
puzzling features of the experiment of Ref. ^ 

We consider a generic system whose Hamiltonian H 
is made of three terms separately describing the multi- 
orbital interacting quantum dot, -ffdot, the conduction 
electrons in the leads, i?icad, and the hopping term be- 
tween dot and leads, H^w: 



H = H, 



dot 



lead 



(2) 



The conduction electrons, in typical experimental setups, 
can freely move in a two-dimensional heterostructure 
[isj . according to the Hamiltonian i?ioad — '^ka ^k'nka, 
where rika — c^.a'^ka and c^^ creates an electron into the 
Bloch state of crystal momentum fc, spin a, and energy 
Ek- The QD Hamiltonian i?dot includes both the single- 
particle and the interacting part: 



dot 



E 



Oil ' ^aicr 



(3) 



In order to label the QD orbitals, we here introduce two 
indices, a — 1 ... , iVdass and i = 1, . . . , iVa, respectively. 
The former index, a, labels the iVdass irreducible repre- 
sentations of the QD point-symmetry group, while the 
latter, i, enumerates the truncated set of Na orbitals 
considered, belonging to the same a representation. We 
assume that QD orbitals form also a basis for represent- 
ing the symmetry of the whole dot -I- leads system [l3| . 
Typical symmetry groups of realistic devices range from 
Dooh to C2V, going from circular to elliptic 0| dots, 
respectively [l^. In Eq. © riaia = clia^aia, c^aia creates 
an electron with spin a in the orbital (a, 1) of energy £0,^, 
and i/int includes the full intra-dot Coulomb interaction. 
Note that Hua does not commute with n^ia except if 
there is either one level only, iVdass — 1, — 1 (non 
degenerate Anderson model), or the Coulomb interaction 
takes an oversimplified form. Here the full inclusion of all 
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Coulomb matrix elements in ffjnt turns out to be crucial 
in the following. Finally, the tunneling term ffmix allows 
for electron hopping between delocalized Bloch states k 
and confined orbitals (a, i) via the matrix elements Vkaf. 



(4) 



kaia 



As a first step, we introduce the zero-temperature ex- 
act retarded Green's function 

iGxx'it) - m (cxW4,(0) + 4,(0)cxW) , (5) 

where either X = aia or X = ka, and (. . .) is the average 
on the interacting ground state |^o) of the whole dot -I- 
leads system in the Heisenberg representation. We work 
with the analytic continuation of the Fourier transform 
of Q in the complex plane of the energy, Qxx'{z)- 

In the following we neglect spin-flip scattering pro- 
cesses [l^ . and therefore we may drop spin indices of Q. 
Moreover, off-diagonal Green's functions of type Gaii3j{z) 
must vanish due to symmetry [but not Gaiajiz), hence- 
forth indicated as G"j{z)]. We first consider the non- 
interacting case (then Hmt = and G is printed in italic) , 
where explicit solutions are readily obtained using the 
Green's function's equations of motions [l3.ll7|: 



fi ^Gkk{z) = 



E 



{z - Ek f 



GUz), (6) 



and 



.)5,m-AZ,Az)]G:,^{z) = h5, 



(7) 



with Ag(z) being the self-energy due to the dot- 
lead interaction: Ag.(z) = J2k^aikVkajiz - £fc)~^ 
In the non-interacting case the diagonal QD Green's 
function assumes the familiar form h^^Gf^{z) = 
[z — Sai — l^fi{z)\~^ , where Im [A^(2:)] is the virtual level 
width, and Re[A^(z)] renormalizes the single-particle 
level £ai- Equation Q is modified to take into account 
electron-electron interaction by introducing the intra-dot 
proper self-energy matrix I]^ (z): 



'Oil ) ^im 



^?„^iz)-AZ,,iz)]g:^^iz)^h6, 



(8) 

Henceforth we focus on the fully interacting system, 
and we consider the case of elastic scattering when only 
a single QD level, (a,T), is coupled to the leads. Indeed, 
this is a reasonable scenario for the electrostatic potential 
barriers separating dot and leads in many experimental 
setups p], where matrix elements Vkai strongly depend 
on both energies Sk and Sai, respectively. We then set 
Vkai = if (a, i) ^ (a,!) and Vfeai ^ with Esa - £cut < 
^k < £ai + £cut , where Ecut is a suitable cutoff. Note that 



Coulomb correlation is included in full and the intra-dot 
self-energy has off-diagonal matrix elements, Sij(z) ^ 
0, a — 1, . . . , iVciass, — 1, ■ • ■ , ^a- This scenario is 
generic enough to correctly describe many experimental 
situations, except the case of degeneracies between QD 
levels £a\ and e^j. 

We proceed in close analogy with Langreth and 
calculate the charge A/ai (in units of e) displaced at the 
dot level (a,!) as the difference between the charge at 
equilibrum in the presence and in the absence of the dot, 
respectively: 



Ma 



1 

TvK 



Im r^"" dz[\Y,Qkk{z) + Q^X^) 



^ E '^kk'^^)^ 



(9) 



where 77 is a positive infinitesimal quantity, /i is the 
equilibrium chemical potential T^, and h^^G^^°{z) = 
{z — ek)~^ is the propagator of a free traveling wave in 
the absence of the QD. By using Eqs. (O and © where 
G's are replaced with C?'s, one obtains the following ex- 
pression for TVqt: 



dz 



1 



9Ag(z) 



dz 



GSiz)- (10) 



We now use the identity 
Im 



dz^^gS{z)^o, (11) 

-r^i (^Z 



which is a natural generalization of the Luttinger rela- 
tion of Fermi liquids '19*1, and it has been already ap- 
plied by Langreth l^j to the single-level case. Combining 
Eqs. ((nil and gives 



TVqt = — Im 



dz—\ng§iz) 
oz 



(12) 



Since the asymptotic form of the Green's function as z — *■ 
—00 is g°j{z) « hSij/z, Eq. (|12ll may be casted into the 
form 



Ms 



1 



Im[lna#(A* + if?)-H 



(13) 



where we have chosen the cut along the positive real axis. 
To connect the displaced charge Aar to the Friedel phase, 
we observe that the T scattering matrix of the single 
channel k is Tkkiz) — \Vkm'^ h~^Q^{z) (cf. Sec. 5.2 of 
Ref. therefore the phase on the energy shell 5F{y,/h) 
is 

Spifi/h) = arg Tfcfc(/i + IT]) = lmlng^{fi + i-q) - tt, (14) 

where we pick the branch of the logarithm as before and 
add the reference constant — tt, so that Q < 5^ < it. 
Comparison of ((T^ and lfn|) gives the desired result: 



SF{fi/h) = ttMs 



(15) 
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The exact sum rule ()15|l states that the Friedel phase 
shift, due to scattering by the dot, is proportional to the 
net charge accumulated at the dot level {a,T) with re- 
spect to the free system (i.e. without dot). In the case of 
a single level (iVdass = 1, -^q = 1), Eq- l(T31) is equivalent 
to the result of Langreth for the Anderson model [l^ . 
In the many-level case, H15|) is a non trivial generaliza- 
tion of previous theories (see e.g. 0| and 0), which can 
be summarized as follows: The sum of phase shifts due 
to consecutive filling of many levels is fixed by the total 
charge which can be placed into the dot, ruled by fj,. If the 
dot levels, possibly broadened due to hybridization with 
continuum states, lie well below fi, and T > Tk, where 
Tk is the Kondo temperature, then the afore mentioned 
charge is an integer quantity fixed by orbital degener- 
acy. According to (|15|l . however, this simple picture is 
generally incorrect, as the following conceptual tunnel- 
ing experiment illustrates. 

Think of varying continuously the chemical poten- 
tial /X across an energy window centered around the 
resonant value hu^cs, which is implicitely given by the 
real part of a certain pole of G^{z), hiOy-cs — ear -I- 
Re [E^(^res) + A§(Sa;res)]- For the sake of clarity we 
here focus on the Coulomb blockade regime only, despite 
the fact that sum rule H15|l applies to the Kondo regime 
as well. The level at fiiUrcs is located between two block- 
aded regions with N and -I- 1 electrons in the dot, 
respectively. The transport window is chosen to be large 
enough to fully contain the width of the QD level, given 
by ^cutj but narrower than the spacing between neigh- 
boring resonant levels. Therefore, as fi is swept upward 
across the QD level, the dot is charged by exactly one 
electron, in addition to those already localized. Nev- 
ertheless, the variation of the spectral density AAai of 
Eq. p5|l across the energy window is generally less than 
one, and consequently the phase shift increment of the 
outgoing wave spreading out from the dot at the top of 
the energy window is less than tt. 

To understand it, note that the total number of scat- 
tering states of the whole dot -I- leads system, in the 
energy window considered above, must be exactly equal 
to the sum of both free travelling waves in the leads and 
confined states in the dots, when the two subsystems are 
decoupled (-ffmix = 0). Therefore, in this case AA/qt niay 
be calculated by simply integrating the spectral density 
of the isolated dot: 



AA/'ai = -— Im 
nn 



t+iv 



dzg\ 



K free 



(^) 



*0 Fail^o 



-£cut+ii) 



(16) 



where |^'^) is the exact interacting ground state of the 
isolated dot with A^ electrons. The following is clear: (i) 
AAai is a positive quantity, which can considerably devi- 
ate from the unit charge that could fill in the level (a,T), 
due to the correlation between electrons localized in the 



dot. AA/ai can even be zero due to spin blockade, namely 
the difference between total spins of j^'o^'''^) and 
is not equal to ±1/2. (ii) The missing weight is recov- 
ered by integration on the whole spectrum, (iii) Equation 
(|16|l provides the basis for exact numerical evalutation of 
the phase shift, by computing the spectral density of a 
confined system e.g . by means of the full configuration 
interaction method |ll[ l2(l l2lj . 

The discrepancy between total and spectral charge 
density of Eq. ((T^ can be neglected only if intra-dot cor- 
relation effects are absent, namely the ground state of the 
isolated dot is a single Slater determinant. In such a case, 
the matrix element appearing in Eq. (|16|l is either one or 
zero, depending if level (a,!) is filled in or not, respec- 
tively, in the A^ ^ A^ + 1 tunneling event. This occurs 
in the degenerate Anderson model [2, ll^ IitI |2^ , where 
H\nt assumes a Hartree-Fock-like form diagonal in the 
(a, i) indices. While such mean-field model is satisfying 
for magnetic impurities in a bulk metal or many-electron 
dots, it breaks down for large dots with very few elec- 
trons, like those of Ref. 0], where correlation effects may 
dominate 0, driving the ground state of the isolated 
dot even towards the Wigner crystallization regime [lll |. 

The idea behind l(T^ is general, as we prove in the 
following case of arbitrary coupling between leads and 
QD levels {Vkai V fc,a,i). We start from the exact 
relation between Friedel phase and delay time TA^y of a 
traveling wave packet when trapped in the dot |23| : 



^ ; = Tdclayl^^j, 



(17) 



where hw is the average energy of the wave packet. We 
assume that Eq. (|17|l . proved for independent particles 
,24] , holds even in the presence of correlation since H^t ^ 
only in the region occupied by the dot, while in the 
outer space the packet moves freely. The delay Tdciay 
can be calculated as the difference between the ( "dwell" ) 
times spent by the packet in f2 in the presence and in the 
absence of the dot, respectively [23ll25|. 

By extending the approach of lannaccone |2(il | , we write 
the wave packet |$) as a superposition of the interact- 
ing incoming scattering states l^'fe), |$} = X]fe (^{^) \^k), 
where \a{k)\^ = 1 so that [$) is normalized to unity. 
The exact stationary interacting state, l^'fc), satisfies 
H\^k) = {Eq + £k) |*fc), with Eo being the energy of 
the ground state without the extra electron to be scat- 
tered, and H |5'o) ~ Eq |*o)- The probability amplitude 
for finding the scattered electron at position r is obtained 
by projectin g 1$ ) on an eigenstate of position of the ex- 
tra electron [20|, vE''^(r,i) |4'o), where we introduce the 
field operator \E'^(r,i), creating an electron in r. The 
field 5'(r,t) may be decomposed onto a mixed basis of 
Bloch free waves (pkif) and confined QD orbitals 0Qi(r): 
*(r, t) = J2k Mr) Cfe + Eai 4>Ur) c^i 113. We therefore 
define the mean dwell time in Q associated with the wave 
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packet |<I>) as 



dt / dr|(*o|«'(r,i)|$) 



(18) 



The integral ^ converges since the probabihty of find- 
ing the scattered electron in Q vanishes for time ap- 
proaching ±00. After decomposing the wave packet on 
the basis made of |^'fc)'s and performing the time inte- 
gration in (|18|l . the delay time is obtained as a sum of 
stationary state contributions |2^, that are separately 
written as 

Tdclay(w) = 2Tr h 5{hlLl - Ek) X 



dr 



|(*o|«'(r)|«'fe)|'-|0fc(r)|' .(19) 



We extend the range of integration in H19|l to the whole 
space, since the contribution outside O is null, and by 
combining (|17|l . (|19|l . and the orbital spectral represen- 
tation of 'l'(r), we obtain the desired result 



1 ddpiu)) 



(20) 



h duj 

where N{lo) is the total spectral density [cf. Q]: 

Oil 

+ ['3^^^'' + i'') - G'^u{^ + i*?)] }■ (21) 

k 

Equation (|20|l is the natural generalization of (|15|l . 
namely a generalized FSR where the phase shift varia- 
tion as a function of energy is proportional to the varia- 
tion of the total displaced spectral density TV. Again, all 
considerations of the previous example [Eq. (fTB|l ] apply, 
i.e. AJ\f does not need to be one between two consecu- 
tive Coulomb blockade regions. The result 1)20(1 is generic 
and independent of the nature of the coupling between 
dot and leads. 

We are now able to focus on the experiment of Ref. |3 
The transmission phase shift of an electron tunneling into 
a QD is determined together with the integer charge 
filling in the dot, which is in the Coulomb blockade 
regime. It turns out that the phase increment between 
specific neighboring conductance valleys, as a function 
of the plunger voltage, is a fraction of tt (■^ between 
iV = 3 and N = A and ~ 37r/4 between N = & and 
N = 7 blockaded regions). Moreover, for these specific 
voltage ranges, the phase variation is continuous and 
smooth. Under such circumstances, 6f and the trans- 
mission phases are expected to coincide IgI ITi . According 
to previous theories based on the FSR |a, Q , one would 
predict a unit increment of tt in all cases, unless some 
additional charge unpredictably accumulates outside the 
dot [2^. Our theory suggests an alternative natural ex- 
planation, i.e. fractional (or even zero) increments of Sf 



originate from strong electron correlation (or spin block- 
ade). This interpretation is supported by the estimated 
low value of electron density. In fact, as the density di- 
minishes, electrons in the dot are expected to crystallize 
|ll| |. which affects dp in a density-dependent manner. Ex- 
tensive numerical simulations will be reported elsewhere. 

In conclusion, we derived an exact generalized Friedel 
sum rule for an interacting multi-level nano-object. The 
variation of the Friedel phase through Coulomb blockade 
regions for values which are fractions of tt is the finger- 
print of electron correlation. 
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Bertoni, and G. Goldoni for stimulating discussions. This 
work is supported by MIUR-FIRB RBIN04EY74. 



[1 
[2: 

[s: 
[4: 

[6; 

[7: 

[8 
[9 
[10: 

[11 

[12 
[13 



[14 

[is: 
[16: 



[17: 
[is: 



[19 
[20 
[21 



Electronic address: lrontani@unimore"itl 

URL: www.nanoscience .unimo . it /max_ index .html 

G. Grosso and G. Pastori Parravicini, Sohd state pfiysics 
(Academic Press, San Diego, 2000). 

A. C. Hewson, Tfie Kondo problem to fieavy fermions 
(Cambridge University Press, Cambridge, 1993). 

G. Hackenbroich, Phys. Rep. 343, 463 (2001). 

S. Datta and W. Tian, Phys. Rev. B 55, 1914 (1997). 
J. S. Langer and V. Ambegaokar, Phys. Rev. 121, 1090 
(1961). 

H. -W. Lee, Phys. Rev. Lett. 82, 2358 (1999). 

T. Taniguchi and M. Biittiker, Phys. Rev. B 60, 13814 
(1999). 

R. Schuster, E. Buks, M. Heiblum, D. Mahalu, V. Uman- 
sky, and H. Shtrikman, Nature 385, 417 (1997). 
M. Avinum-Kahsh, M. Heiblum, O. Zarchin, D. Mahalu, 
and V. Umansky, Nature 436, 529 (2005). 

C. P. Garcia, V. Pellegrini, A. Pinczuk, M. Rontani, 
G. Goldoni, E. Molinari, B. S. Dennis, L. N. Pfeiffer, 
and K. W. West, Phys. Rev. Lett. 95, 266806 (2005). 
M. Rontani, C. Cavazzoni, D. Bellucci, and G. Goldoni, 
J. Chem. Phys. 124, 124102 (2006). 

D. C. Langreth, Phys. Rev. 150, 516 (1966). 

The specialization of the formalism to two separate leads, 
as it is usually done in the QD literature, is straightfor- 
ward but irrelevant in this context. 

A. P. Klein and A. J. Heeger, Phys. Rev. 144, 458 (1966). 
G. Bester and A. Zunger, Phys. Rev. B 71, 045318 
(2005). 

Since most of the point-symmetry groups of realistic QDs 
are abelian, the Na orbital levels belonging to the same 
irreducible representation a are generically non degener- 
ate. 

P. W. Anderson, Phys. Rev. 124, 41 (1961). 
A. A. Abrikosov, L. P. Gorkov, and L E. Dzyaloshin- 
ski. Methods of quantum field theory in statistical physics 
(Dover, New York, 1975). 

J. M. Luttinger and J. C. Ward, Phys. Rev. 118, 1417 
(1960). 

M. Rontani and E. Mohnari, Phys. Rev. B 71, 233106 

(2005) . 

M. Rontani and E. Mohnari, Jpn. J. Appl. Phys. 45, 1966 

(2006) . 



5 



[22] A. Levy Ycyati, F. Floros, and A. Martm-Rodero, 

Phys. Rev. Lett. 83, 600 (1999). 
[23] F. T. Smith, Phys. Rev. 118, 349 (1960). 
[24] E. P. Wigner, Phys. Rev. 98, 145 (1955). 
[25] E. H. Hauge and J. A. St0vneng, Rev. Mod. Phys. 61, 

917 (1989). 

[26] G. lannaccone, Phys. Rev. B 51, 4727 (1995). 



[27] We here ignore any complication due to the non or- 
thogonality of (j>k{f) and <f>ai{r). For a discussion see 
e.g. T. E. Feuchtwang, Phys. Rev. B 10, 4121 (1974). 

[28] A. Levy Yeyati and M. Biittiker, Phys. Rev. B 52, 14360 
(1995). 



